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We study the entanglement properties of random pure stabilizer states in spin-^ parti- 
cles. We obtain a compact and exact expression for the probability distribution of the 
entanglement values across any bipartite cut. This allows for exact derivations of the 
average entanglement and the degree of concentration of measure around this average. 
We also give simple bounds on these quantities. We find that for large systems the 
average entanglement is near maximal and the measure is concentrated around it. 
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1 Introduction 

Entanglement is a fundamental resource for quantum information processing. The classifi- 
cation and quantification of this resource in two- and many-partite systems is therefore of 
significant concern, as can be seen in the review papers [U [2l [S] IH [5l [6l [7] . Over the past 10 
years or so the properties of bi-partite entanglement have been explored in some detail and 
many of its basic features are now reasonably well understood. However, the entanglement 
properties of multipartite systems are far more complex and our current understanding of this 
setting is limited. 

There are various approaches that one might take to achieve progress. Firstly, one may 
impose additional constraints on the set of states and/or the set of operations that one is 
interested in, ideally without significantly reducing the variety of possible qualitative entan- 
glement structures. In this context an interesting class of states that arises is that of stabilizer 
states. Despite having various restrictions, these stabilizer states possess a rich entanglement 
structure exhibiting multi-partite entanglement [8l [9l fTO l [TT | [T2 \ [T3] . 

With increasing numbers of particles the finest classification of entanglement results in 
a number of types of entanglement that grows exponentially with the number of particles. 
However, many of these types may be atypical in that their probability, with respect to some 
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natural measure on the set of states, vanishes in the hmit of large numbers of particles. 
This suggests a different approach that studies only the entanglement properties of typical 
states. Even before the emergence of quantum information such questions had already been 
of interest. An example is the study of the expected entropy [Ml [HI [16] of a subsystem when 
averaged over the invariant measure of pure states. In quantum information this entropy has 
received considerable recent attention as it is precisely the entanglement between this subset 
and the remainder of the system. The mean entanglement has been studied, as well as the 
probability for deviations from this mean. These can be shown to be exponentially decreasing 
with the difference from the mean, a property known as 'concentration of measure' [18j . 

In the present paper we are interested in a combination of these two approaches, i.e. we 
study the typical and average entanglement of random stabilizer states. These depend on 
the probability distribution of entanglement, which gives the probability of finding a certain 
entanglement in a stabilizer state that has been picked at random. Therefore the objective 
here is to find and study the entanglement probability distribution of randomised stabilizer 
states. 

The paper is organized as follows. In section 2 we provide some basic results concerning 
random states, define stabilizer states, the basic technical tools and results that will be used 
subsequently. In section 3 we state and prove our main result. We study the probability 
distribution of the entanglement across any given bipartite split of the system. We provide 
a compact, explicit and exact formula for this probability distribution and present its proof. 
In section 4 we use this result to study the average entanglement of a set of spins versus the 
rest and demonstrate that this distribution implies a concentration of measure, i.e. for large 
numbers of particles the probability that the entanglement of a specific state will deviate from 
the mean value decreases exponentially with that deviation. 

2 Basic techniques and definitions 

In the following we present some basic results concerning general random quantum states as 
well as basic tools for the description of entanglement in stabilizer states. 

2. 1 Entanglement of randomised quantum systems 

Consider a system of N spin- 1/2 particles and pick random quantum states from the unitarily 
invariant distribution. This may for example concern, in an idealized setting, a gas of two-level 
atoms. As the atoms collide and interact at random their energy levels become entangled. 
Asymptotically the distribution on pure states becomes uniform, such that any pure state is 
equally likely. For such a distribution of states we may ask: What is the average entropy of 
entanglement WjSa{Nai Nb) of a set of Na spins? Some of the first studies of this question 
were [T3J[T7] and the explicit solution('Page's conjecture') was conjectured in [13] and proven 
in [16] . The explicit solution is given by 



with the convention that Na < Nb and where Na + Nb = N, the total number of particles. 

This can be used to show that the average entanglement is very nearly maximal, meaning 
close to Na, for large quantum systems, i.e. iV ^ 1. Hence one concludes that a randomly 
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chosen state will be nearly maximally entangled with a large probability. Indeed, it was 
recently shown that the probability that a randomly chosen state will have an entanglement 
E that deviates by more than 5 from the mean value WjSa{Na, Ng) decreases exponentially 
with 5^. This phenomenon, alongside many other properties, which is known as concentration 
of measure of entanglement, was proven and extensively studied in [T5] . 

2.2 Entanglement of stabilizer states 

As it is our aim to study the typical properties of stabilizer states we use this subsection to 
present a number of basic tools and observations that are useful in this context. 

Stabilizer states are a discrete subset of general quantum states, which can be described by 
a number of parameters scaling polynomially with the number of qubits in the state [H [9l I22j . 

A stabilizer operator on N qubits is a tensor product of operators taken from the set of 
Pauli operators 

and the identity /. An example for = 3 would be the operator g = X ^ I ig) Z . A set 
G — {<?!, ■ • ■ ,gK} of K mutually commuting stabiliser operators that are independent, i.e. 
Y\f=i ~ ^ exactly if all Si are even, is called a generator set. For K — N a generator 
set G uniquely determines a single state {tp) that satisfies 5fe|V') = IV') for all /c = 1, . . . , A^. 
Such a generating set generates the stabilizer group. Each unique such group in turn defines 
a unique stabilizer state. 

A first observation that will be useful for the following considerations is the fact that the 
bipartite entanglement of a stabilizer state, i.e. the entanglement across any bipartite split, 
takes only integer values [TOl [12] . 

In the proof of our main theorem in section III below we will furthermore use results from 
[12]. In particular we use their result 1, that a stabilizer state of E ebits can be generated 
from: 

• 'local generators' of type gA ® Ib and Ia® gB where g refers to a member of the Pauli 
group. 

• 'non-local generators' of type gA®gB- These generators come in pairs where the entries 
corresponding to Alice (and Bob respectively) anti-commute. 

The subgroups generated by the local elements are labelled Sa and Sb respectively, and the 
subgroup generated by the non-local elements Sab ■ The entanglement E oi a. state is given 
by the number of pairs in the minimal generator set of all non-local pairs, i.e. 

E=\Sab\I2. (3) 

For example the GHZ state |000) -I- |111) with respect to the division Na = 2 (first two 
particles) and Nb = 1 is defined by the generator set {XXX, IZZ, ZZI). The local generators 
on Alice's and Bob's part respectively are given by {ZZI) and (/). The non-local part is 
{XXX, IZZ) which consists of only one pair so that the entanglement \s E = 1. We will also 
use Eq. 6 from ^12j giving entanglement E as 



E = Na-\Sa\=Nb~\Sb\ 



(4) 
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where the | . | refer to the size of the minimal generator set of the subgroup in question. 

FinaUy we use the fact that in order for the stabiHzer state to be non-trivial it is necessary 
and sufficient that the elements of the stabilizer group (a) commute, and (b) are not equal to 

-im- 

3 Main Result- Entanglement Probability Distribution 

Let us consider the setting shown in fig. 1 where a system of N spin-^ particles is split into two 
arbitrarily sized subsets consisting of Na and Nb particles such that N = Na + Nb ■ Consider 
now a randomly chosen pure stabilizer state on the N particles which is chosen according to 
the uniform measure. Generally the spins in set A will exhibit entanglement with the spins 

ooooooooooo o 

Fig. 1. The spins are grouped into two sets A and B. Set A contains Na spins, set B contains 
Ng and the total number of spins is iV = Na + Ng . 



in set B and this entanglement will take discrete values e {0, 1, 2, . . . , min{NA, Nb)}- The 
main result of this paper is an exact and compact expression for the probability distribution 
of the entanglement values E under randomised stabilizer states. 

Theorem I —Entanglement probability distribution: In a system of N spins where Na{Nb) 
is Alice's (Bob's) number of qubits the probability of finding that entanglement between A 
and B equals E in a randomly chosen pure stabilizer state is given by 

_ UlM^' + 1) A (2^-^^+^-^- - 1) (2^^+^- - 2^^--^) 

^ ' (2^ + 11-^^ 22J -1 

where E is an integer and < E < min{NA, Nb)- 

Proof: The strategy of the proof is to firstly define the probability distribution on stabilizer 
states. Then we count the number of states for which the entanglement between sets A and 
B is given by E and label this number by n^. Finally the probability weight is multiplied by 

to obtain the probability distribution of entanglement. 

For a set of N spins we denote the total number of possible stabilizer states by ntot{N). 
Randomising these states is then achieved by the random application of unitary maps that 
take any stabilizer states to another stabilizer state. The probability distribution that is 
invariant under randomisation is the uniform distribution. Therefore the probability of picking 
a stabilizer state j at random is given by pj = 1{n) ■ 

For given values of N , Na and Nb, let E denote the amount of entanglement between set 
A and B, and let n^^^\N) denote the number of stabilizer states on N particles realizing 
this value of entanglement. It will be our task to determine n^^^\N) exactly. It is clear that 
for any Na it must be the case that ntot{N) satisfies: 

lN/2} 
_B=0 
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where \x\ denotes the largest integer smaUer than x. 

The proof of our theorem wiU now be presented in the form of several useful lemmas. 

Lemma I —Total number of states: The total number of distinct stabilizer states for N 
particles is given by 

N 

ni„t(iV) = 2^ J](2' + l)- (6) 

fe=i 

Proof: This is proven in [llj and an alternative proof employing different techniques is found 
in [2n]n 

Lemma II —Number of separable states: The number of pure stabilizer states for which 
the spins of sets A and B are not entangled is given by 

n[''^\N) = ntot{NA)ntot{NB). (7) 

Proof: A pure state does not exhibit entanglement between the spins of sets A and B exactly 
if it is a tensor product between the two sets. Therefore the total number of unentangled 
pure stabilizer states is the product of the total number of states on each partD 

Lemma III —Invariant ratio of probability distribution: For < E < min{NAT Nb) and 
N ^ Na + Nb we have 



n 



E V ^ -^J _ "-E 



(8) 



Proof: This is true exactly if there is, for arbitrary Na and < -E < min(NA, Nb), a 
constant ratio 

n(f-+^)(iV+l) 



n^E^\N) 



f{E,NA) (9) 



'E 

where f{E,NA) is some function that only depends on E and Na but not on Nb- To see 
this we now employ the methods described in section II. Note firstly that the numbers he are 
defined by counting how many stabilizer states there are corresponding to a given entangle- 
ment. This is given by the number of distinct stabilizer groups with that entanglement. For 
this number to be non-zero we only consider counting it when < E < min(NA, Nb)- The 
stabilizer groups have to be Abelian and contain no negative identity, as mentioned in the 
introduction. Consider now the number of ways of picking stabilizer group elements, under 
these restrictions, on Bob's side. We can call this number /3 and the corresponding number on 
Alice's side a. The problem is simplified by counting the number of distinct generator sets of 
the group. Several generator sets will generate a given group, so the number of generator sets 
is larger than the number of groups by some overcount factor. We now show that (3, as well 
as the overcount factor, are invariant during the transition r : A^^ i—s- Na + 1 under constant 
entanglement E. 

Firstly note that the contribution of Bob's side to /3, and to the overcount factor, is 
independent of \Sa\- This is because there are only identities in the corresponding terms on 
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Bob's side. It does depend on \Sb \ and \Sab\ though. Eq. [4]iniphes that \Sa\, but not \Sb\ 
changes under the transition t : Na ^ Na + 1- The condition of constant E is equivalent to 
constant \Sab\ by Eq. [3l 

So r : /9 I— > /3 and t : a i-^ a' . Therefore Bob's contribution to the ratio ([9]) cancels; the 
ratio is independent of Nb- □ 

Now we note that the following closed form expressions 

oNa+E 92B-1 , 

n(f-)(Ar) = ^2^+^-Na-e 1) nZiiN) (10) 

satisfy equations ([5]), ([6]), ([7|) and ([8|), for < £' < min{NA, Nb), i.e. for all possible settings. 
What remains to be shown is that this choice is indeed unique, i.e. that we can construct 
'n}j^^\N) recursively and uniquely from equations (O, ([S]), ([7]) and We proceed in the 
following steps. 

• Note that we know n^j^^\N) for Na = 1 and all valid choices of for = by virtue 
of Lemma II and then for _B = 1 by virtue of eq. ([5]) . 

• Assume that we have obtained n^^"^\N) for Na — 1, 2, . . . ,ro and all valid choices of 
E € {0, . . . , Na} and N € {Na + E, . . .}. Now we will demonstrate that this uniquely 
defines np^'^^N) for all £; e {0, 1, . . . , ro + 1} and all N e {ro + 1 + E, . . .}. 

• To this end realize first that n^E°X^\-^) again known for all N by virtue of 
Lemmas I and II. Now assume that for fixed ro + 1 we have found ni''°^'^\N) for 
X G{0,1,...,E} and all N. 

• We find n'"^'^^\N) for all in the following way. First we use eq.® to obtain 

n^ElV\E+ro + 2) = ntot{E+ro + 2)-J2f^nn^P^^\E+ro + 2). Now we employ the 

recursion relation in Lemma III to obtain n'^'Tt"'^'' (A^+1) — "^/t\ np^l (N) 

where each term on the right hand side is already known by assumption. 

This completes the construction for arbitrary TV which in turn completes the construc- 
tion for arbitrary E. Therefore, all n^^'^^N) are uniquely determined. It is now cum- 
bersome but straightforward to check that eg. pO]) satisfies all the recursion relations 
which in turn implies that it is the unique solution presenting the correct values for 
n(f-)(iV). 

This finishes the proof for Theorem I □ 
4 Implications 

This section will present a discussion of the implications of Theorem I. 
4-1 Product-free form 

Theorem I can be more conveniently evaluated in the format free from products 
Corollary I —Product-free form: The probability of E entanglement is 
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where Na, Nb and N are as defined earlier, and 




Si 



(12) 




and finally 



E 



S2 



) + log^il - 2-^+^--i+^-) - log^il - 2-^') . 



(13) 



fc=i 



Proof: This is a reexpression of Theorem I □ 

We will see Si and S2 converge to numbers of magnitude <C N when N is large so they can 
be seen as minor modulations to the leading Gaussian-type behaviour. This statement will 
be made more rigorous later on in this section. 

4-2 Comments on entanglement probability distribution 

It is instructive to consider qualitatively what kind of distributions Theorem I describes, before 
making precise mathematical arguments. Firstly consider the simplest non-trivial case, that 
of N=2, Na = Nb = 1. Here one obtains that the number of states with one ebit is ni = 24, 
whereas uq — 36 with ntot — 60. Hence the expected entanglement is 0.4. If we now increase 
Bob's size to Nb — 2, the formulae show that tiq = 360, ni = 720 with ntot = 1080, giving 
an average entanglement of — |. When considering all possible Na under increasing N, 
one observes that the maximal and near maximal values of entanglement become the main 
contributors to the number of possible states. This can be seen from Eq. (fTTj) in which the 
leading terms inside the sum correspond to a Gaussian centred on N/2. Note however that 
although the corrections to this Gaussian, such as the E2 term are small, they contain an 
_E-dependence that adds some subtleties to this behaviour that will require a more rigorous 
analysis. This can be understood by noting from Eq. (fTO)) . that 



where the second term quickly disappears with increasing N. If Na is constant, the ratio 
tends to exponentially fast. However if Na — y then UN^ -^/riNji 2 for large N which 
implies in fact that the maximum of the probability distribution is shifted to N/2 — 1. In that 
case one sees from Eq. (jlOp that n^j^ ^/rif^^ -^ — > |, after which the ratios are very small. 
Hence for this case Na — 1 is the most likely entanglement of a state picked at random rather 
than Na- 

We can now use the entanglement probability distribution to derive the average entangle- 
ment of randomised stabilizer states. 

4.3 Average entanglement of randomised stabilizer states 

Corollary II —Average entanglement: The average entanglement, ES"^, in stabilizer states 
sampled at random is given by 



1 1 



(14) 



2N-2Na _ 1 2^ — 2^2iVA-i 
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where Na is the number of qubits belonging to Ahce and N is the total number of qubits. 
We follow, without loss of generality, the convention that Na < N — Na = Nb- The total 
state is bipartite and pure and the average is taken over a flat distribution on stabilizers. 
Proof: We take as a starting point 

Na 

E,Sa = P{E)E (15) 

E=l 

which is then evaluated using Theorem IQ 

It is also worth mentioning that the case Na — 1 yields the simple form 

ntot-na 3 

4-4 Concentration of measure 

If some property becomes very likely in certain circumstances, we say the measure is becoming 
concentrated around states with this property. It is known that the measure on general 
pure states is concentrated around states with near maximal entanglement. Here we show 
that a similar concentration takes place on stabilizer states, with some subtleties. We know 
that the mean is close to maximal which we can use as an upper bound. Below we give an 
exact expression for the probability of picking a state with an entanglement less than the 
average by some number. We then give an inequality for the same quantity. 

Corollary III —Concentration of measure around mean in stabilizer states: The proba- 
bility of picking a stabilizer state such that it is at least e less than the average ESa is given 
by 

P{SA<TESA-e)^ PiE) (17) 

E=l 

where P(E) is given by Theorem I and WiSa by Corollary II. e is a constant which can take 
any value in the range of 1 to Na — 1. 

The concentration of measure can also be given in the shape of an inequality by exploiting 
the fact that the leading behaviour of the probability distribution of entanglement is Gaussian: 

J2 P{E)<i{N,Na) di;e-'"2W2-ijf (^8) 

E=i 



where 



7(iV, Na) = 2 ^ ^ ' ■ (19) 



si is defined in Eq. ([22|l . S2 is given by 

hNA t uNb 

S2 = ——, ;o92(15/16) + /i (20) 

where for typographical reasons — /o(72(3/4) — h. 

Proof: The terms inside the integral are taken from the product free form of the prob- 
ability distribution. We then factor out 7. This is possible by using i?-independent bounds 
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on Si and E2. The latter is an £'-independent version of Eq. [26] below. The integral can be 
taken as a bound to the sum, as the remaining Gaussian is centred on A^/2, and without loss 
of generality Na < N/2, meaning there is a monotonic rise within the bounds of the integralQ 
The Gaussian integral can be evaluated using the error function. 

4-5 Comparison with general random states 

The typical entanglement of stabilizer states is very similar to that of general states. By 
this statement we refer to two properties of the entanglement in random large general states 
picked from the unitarily invariant distribution: the average entanglement is near maximal, 
Na, and the probability distribution is exponentially concentrated around this average. 

The concentration of measure section above shows that the latter is true also for random 
large stabilizer states. The typical entanglement of stabilizer states is nearly maximal, as is 
the case for general states. The averages are slightly different though, and the concentration 
of measure around the average appears to be a bit less abrupt in the stabilizer state case. 

Figure 2 shows how the averages of randomised general states and stabilizer states compare 
for the full range of possible Alice and Bob divisions where the total system is ten qubits. 



Averages of general and stabilizer entanglement 
are very similar 




Fig. 2. For a system of ten qubits with Na belonging to Alice the average entanglement of both 
general [15] and stabilizer states (Corollary II) are shown. 

In order to prove this similarity for arbitrary Na a good starting point is a rearrangement 
of Corollary II. 

^Sa = Na{1 -^)-y (Na- ij)2^^^^^-(f -^)^+^^+^^ (21) 
ntot 
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which uses the product free form of Theorem I as in Eq. To obtain a lower bound on 

this we will need the following upper bounds on Si and S2. 
•Si upper bound: 

A compact expression for Na — 1 was given in eq. (|16p so that we may concentrate on the 
case Nji > 2. Then Si < si where 

Na-2-\ 

Si =Zo52(3/2) + ?og2(5/4) 



13 
36 



1 



(22) 



Proof: Firstly note 



Na N 

Si = ^ ^032 (1 + 2-^1) - 1092(1 + 2-^') 

]1=1 j3=NB + l 



Na 



< ^°92 (1 + 2- 



?og2(3/2) + ?o.g2(5/4) 



Na 
J 1=3 



(23) 



(24) 



(25) 



Now consider using an upper bound of the type exp{a + bk) on log2{i + 2 ^ ''), where a and 
b are constants to be determined. Choose a and 6 such that the bound is exact for k — I 
and fc = 2 respectively. This yields a = 2ln{log2i9/8)) - ln{log2il7 = -1.10825... and 
b = /n(/og2(17/16)) - ln{log2{9/8)) ^ -0.664143.... That exp{a + bk) is indeed an upper 
bound to log2{l + 2^^^'') for k > 2 follows by immediately comparing the gradients of the 
two functions. One now substitutes this bound into Eq. (j25p and uses the standard formula 
for sums of geometric sequences, as well as bounding the messy-looking logarithm ratios by 
rational numbers, to recover Eq. (|22|) n 

The upper bound on S2 can be derived employing very similar ideas. 
• S2 upper bound: To achieve a more compact notation let h = -^0(72 (3/4) 



S, < 



1 - h 



-E 



^-2E 



1) Zo.g2 (15/16) + 



(26) 



1 - /i J 3 

Proof: Essentially the same method is used as in the bound for Si. We use two inequalities 
that are valid within the range used only. Firstly 

- log2 (1 - 2"2-''"2) < e'"(-'°92(15/16))-2(n2+2in2j ^27) 

which was found using an exp(a + bj) type bound which we required to be exact for j — 1. 
The maximum allowable gradient such that it is an upper bound for j >1 was used. Secondly 

Jn(~log2(Z/A)){l-x) ^28) 



'052 (1 



< 



which was found using an exp(a + bx) type bond, where a and b are chosen such that the 
bound is exact for x — 1 and x — 2. The gradients of the two functions are such that it will 
be an upper bound for all points in the range. These two bounds are then inserted into the 
sum in the definition of S2 and with the standard formula for sums of geometric sequences 
one recovers Eq. ([^ □ 

Using these bounds in Eq. (pT|) provides a very tight fit to the exact result. 
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5 Summary and Outlook 

In this paper we have studied the entanglement properties of randomised stabihzer states on 
N particles. We considered a bi-partite setting with two sets containing Na and Nb spins 
with N = Na + Nb and the entanglement between them. We obtain an exact and compact 
probability distribution for the possible values of entanglement (which are integers) and use 
this to give exact values for the average entanglement and the concentration of measure 
around it with increasing size of the system. It was found that, for large systems, the average 
entanglement is nearly maximal and the probability of a state having a significantly different 
entanglement than this average decreases exponentially with the difference. We also compared 
these results to the case of general random quantum states and found a close similarity. 

This suggests that it would be interesting to conduct research into how far this similarity 
carries, in particular in the context of multi-partite entanglement between large and/or non- 
contiguous sets of particles and of dynamical features of quantum states under the action of 
random quantum gates (stabilizer or general unitary gates) |171 121j and their approach to 
the equilibrium distribution. These topics will be the subject of a forthcoming publication 
currently in preparation. 
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